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Abstract. We argue, and find through numerical work, that the results 
of non-dynamical Monte Carlo computer simulations cannot be applied to 
describe the formation of topological defects when the correlation length 
at the Ginzburg temperature is significantly smaller than the horizon size, 
the case which was originally considered. To test the current hypothesis 
that ‘infinite’ strings at formation are essentially described by Brownian 
walks of size the correlation length at the Ginzburg temperature, we ex- 
amine equihbrated fields at the Ginzburg temperature. We find that no 
‘infinite’ structure exists in equilibrium for reasonable definitions of the 
Ginzburg temperature, and that horizons must be included in a proper 
treatment. A phase transition, from smaU-scale to large-scale string or 
domain wall structure, is found to occur very close to the Ginzburg tem- 
perature, in agreement with recent work. We also investigate the forma- 
tion process of domain walls and global strings through the breaking of 
initially ordered states. To mimic conditions in the early Universe, cool- 
ing times are chosen so that horizons exist in our sample volume when 
topological structure formation occurs. The classical fields are evolved in 
real-time by the numerical solution of Laugevin equations of motion on a 
three dimensional spatial lattice. Our results indicate that it is possible 
for most of the string energy to be in small loops, rather than in long 
strings, at formation. 
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CA 95064. 


Optrated by Universities Research Association Inc. under contract with the United States Department of Energy 


- 1 - 


I. Introduction 

About a decade ago, some important implications of cosmological pha^e transitions 
were realized. The possibility that a discrete symmetry in the early Universe be- 
came broken, resulting in domain structure, was discussed by Zel’dovich et al [1], 
Shortly thereafter. Kibble [2] classified topologically stable defects, and discussed 
their subsequent evolution in the Universe. To review, the essential idea behind 
defect formation is that below a critical temperature degenerate vacuum states are 
realized. The most direct argument for the formation of such defects is then based 
on causality- which requires that different horizon volumes end up in uncorrelated 
vacuum states. Depending upon the symmetry that is being broken, regions of the 
old, symmetric vacuum state can be trapped in point defects (monopoles), line de- 
fects (strings), and planar defects (walls) as a consequence of neighboring horizons 
being in different vacuum states. More complicated objects, such as walls bounded 
by strings, and monopoles connected by strings can also form. In the typical field 
theories describing these objects, their continued existence in the Universe, after the 
phase transition, is purely due to non-thermal effects. That is, the ground state of 
the theory at zero temperature is where all of space is in just one vacuum state, and 
there are no topological defects- an unattainable state, however, as long as there are 
horizons. 

Of all the defects, strings have received special attention as they can produce the 
density fluctuations required for galaxy formation [3], provided the energy per length n 
is GpL w 10~® (a typical GUT scale). The success of the scenario hinges upon whether 
or not a scaling solution exists, i.e., whether or not long strings can efficiently lose 
energy through loop production to avoid a string-dominated universe. Much work 
has gone into numerical simulations of Nambu strings in an expanding Universe [4] to 
determine if such a solution exists-the evidence suggests there is a scaling solution. 
To better understand the properties of strings, the distribution of strings at formation 
has been modeled by using a non-dynamical Monte Carlo simulation-which is also 
used as initial conditions in the simulations. In addition, equilibrium properties of 
strings has been studied via analytic treatments of string statistical mechanics and 
by computer simulation. However, all previous models of formation lack a notion of 
dynamics, an essential ingredient of the formation process. 

The formation of cosmic strings has been extensively studied using non-dynamical 
Monte Carlo simulations [5,6,7] . In these simulations, space is discretized on a lattice 
and the phase x of a complex Higgs field $ (= /e’^, where / and x are real) is randomly 
assigned values [0,27 t] to the lattice sites (or, alternatively, to the cells) which represent 



- 2 - 


regions of constant phase. A frequent approximation, which is really not necessary, is 
to discretize the infinite number of vacuum states into a few possibilities (e.g. x = 0? 
27 t/ 3, 4?r/3). Changes in phase Ax between neighboring sites are taken minimally, as 
gradients in the Higgs field cost kinetic energy. Upon traversing any closed path, the 
net phase change Ax is 2irN, where N is an integer (the winding number). Strings, 
which must either be closed or have ends on the spatial boundary, are then easily 
found by checking each lattice face for a non-zero winding number. On a cubic 
lattice, following the above procedure for the identification of a string, (7V( <1. If 
more than one string enters (and exits) a cell, the incoming and outgoing strings 
are randomly connected. The results of such a simulation indicate that ~ 80% of 
the string length is in long strings (> simulation dimensions), with the remaining 
~ 20% in ‘small’ loops with a scale-invariant distribution. For reference, and further 
comparisons, we plot the distribution of strings obtained using this technique, for one 
run on a 35^ lattice, in Fig. 1. 

For the above model to be self-consistent, the lattice sites must be separated by 
a length A over which the phase x is uncorrelated. The horizon length provides an 
upper bound to A. However, if one takes A to be the horizon length, then the large- 
scale structure of the string distribution is obtained {assuming x does not fluctuate 
too much on scales << A), but all information on sub-horizon structure is lost. It has 
been suggested [5] that the full string distribution, at formation, could be obtained 
by identifying the lattice spacing A with the correlation length determined from 
spatial variations in $ at the Ginzburg temperature Tq- The strings then essentially 
have the shape of Brownian trajectories of persistence length formation [5]. [For 
a good review of the heuristic picture of string formation, see Ref. [8].] However, we 
point out that it is not clear that their interpretation of the lattice spacing as is 
consistent with their model of randomly assigning phases. Troublesome points are: 
(1) does not soley describe correlations in X! if ^■Iso reflects spatial fluctuations 
in the magnitude of $, i.e., / (2) there are two physical length scales that arise 
in this problem, the correlation length ^ and the horizon size d, yet their picture 
only involves ^ (3) since horizons do not explicitly enter in their model, we might 
extrapolate their results to describe equilibrium configurations of strings; however, if 
the fields in a very large volume were completely thermalized, the symmetry should 
be broken at the Ginzburg temperature-in contrast to what is obtained by randomly 
assigning phases with equal probability. The formation picture given by Vachaspati 
and Vilenkin is very appropriate for the case — d, yet this is not the case they 
had in mind (^g << d), nor does it represent an interesting dynamical range in the 
parameters d. It is certainly plausible that the horizon length scale could have 
little effect on the initial string distribution for ^g << d, e.g. if thermal fluctuations 
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in the Higgs phase are sufficient at the Ginzburg temperature to produce large-scale 
structure, but this must be checked (we shall find that the horizon plays an extremely 
important role). 

To elucidate the details of topological structure formation, we numerically solve 
finite-temperature equations of motion, describing both domain walls and global 
strings, with additional noise and damping terms which represent coupling to a heat 
bath. The equilibrium defect distribution is examined at various temperatures and, 
in particular, at the Ginzburg temperature To- In addition, we model dynamics by 
preparing ordered states and then cooling to Tq, where the structure is examined, 
on time scales that allow horizons to be present in the sample volume. We, there- 
fore, simulate, and numerically verify, the ‘Kibble mechanism’ for the production of 
topological defects. A number of questions are also addressed. For example, what 
qualifications, if any, should be placed on the presently held view that 80% of the 
length in strings at formation is in long strings? Is the only relevant scale in de- 
termining the string distribution, as has been suggested? If the horizon does in fact 
play a role, what does the string distribution look like at formation if the correlation 
length at the Ginzburg temperature is much smaller than the horizon? 

The paper is organized as follows: in Sec II we provide the details of our approach 
and examine the properties of domain walls at formation; in Sec III we examine the 
properties of global strings at formation; in Sec IV we summarize our work and make 
some concluding remarks. Both Sec II and Sec III are subdivided into equilibrium (no 
horizons) and non-equilibrium (horizons) sections. The equilibrium sections test the 
standard hypothesis that infinite strings or walls can exist at the Ginzburg tempera- 
ture even when there are no horizons. In the non-equilibrium sections, horizons are 
introduced to investigate their role in the formation process, as we shall determine in 
the equilibrium sections that the standard hypothesis needs modification. 


II. Formation of Domain Walls 

A. Equilibrated Fields 

We first consider the simplest topological defect, the domain wall, which is character- 
ized by one real scalar field a and a zero-temperature bistable potential which 

we take to be V{a) — A(<r^ — t;^)^/ 4!. At finite temperatures, the effective potential 
V’r{<t>)i where <f> is the cleissical part of <r, becomes relevant. The exact nature of 
the temperature-dependence in the model is not very important here, and we sim- 
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ply choose a Ginzburg-Landau type potential with Vj{<f>) = —m^T<f>'^f2 + A(^^/4! + 
constant, where r = 1 — {TjTcY and Tc 2r} (m and A are taken to be constants). 
We could just as easily take a phenomenological viewpoint and consider the form 
originally studied by Landau (r = 1 — TfTc), but choose the other form as it repre- 
sents the leading temperature dependence in the high temperature expansion of the 
full one loop finite-temperature potential calculated from V(<r), from which one finds 
Tc ~ 2r] [9]. [To be consistent with this approximation, that the temperature be sig- 
nificantly larger than the mass, we stay reasonably close to the critical temperature 
in our simulations.] At temperatures above the critical temperature Tc, an ordered 
state < (f> >= 0 exists, but as the temperature is dropped below Tc the <p ~4> 
discrete symmetry will be broken, and at T = 0 either < <}> >= —t) ot < <f> >= +t}. 
In the early Universe, however, thermal contact is limited to a horizon volume. This 
means that the (f> *-* —<t> symmetry will be restored on scales larger than the hori- 
zon. The regions that separate the two distinct vacua are domain walls, and their 
large-scale existence is guaranteed as long as complete thermalization cannot take 
place, e.g. there are horizons. The classical field theory solution of a domain wall in 
the x-y plane at zero temperature is <^ = r/tanh( 2 m/\/ 2 ), which yields a wall width 
w ^ \/2m~^ [m = y/Xrj/y/2, the scalar mass]. Operationally, the classical solution is 
a useless definition of a domain wall. We shall simply define a wall to be a surface 
with (^ = 0, at any temperature. 

To study the formation of domain walls, we solve the (flat space) classical equation 
of motion with two additional semi-phenomenological terms F<^ and 

^ _ vV + dVT{(f>)ld<i> = C (2.1) 

where dots denote time derivatives, and the field ^ is a noise term, controlled by 
an amplitude A, which represents coupling to a heat bath and has the following 
properties; 


< C{x,t) >- 0 

(2.2) 

< C(a:,t)^(x', t') >= 2TA8^{x' — x)8{t' — t) 

(2.3) 


is a dissipative term in the equation of motion, which must be present if finite- 
temperature equilibrium configurations are to be established. In an expanding Uni- 
verse such a damping term arises naturally; however, without loss of the physics at 
hand, we do not consider expansion and we choose F to be a constant. Expansion 
would have several qualitative effects: (1) damping, given by F = l/2t in a radiation 
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dominated Universe, (2) horizons of size 2t, and (3) redshifting of the V^<f> term in 
Eqn. (2.1) by a factor oc Since we are presently testing the standard argument 
that infinite structure can result from the equilibrium distribution at the Ginzburg 
temperature, ignoring horizons, we certainly do not care about the expansion effects 
(2) and (3). The only crucial property in this investigation is damping (so equilibrium 
configurations can be achieved), which need not even be supplied by the expansion 
of the Universe. Damping will also occur due to the frictional effects associated with 
the couphng of <j> to other fields. On dimensional grounds, we simply take F ~ rUg//, 
where the effective mass mg// = m-y/r. The choice of F, for a fixed temperature, will 
only affect the probability distribution of the (j> field when we consider non-equilibrium 
field configurations (see below). In any case, the damping rate chosen here will not 
significantly differ from the expansion rate for parameters of interest. 

Eqn (2.1) is analytically insoluble, even without the additional terms, so we must 
numerically solve a discretized version of the equation of motion. We obtain the 
equations of motion from the lattice Hamiltonian H, which we take to be: 

+ A V i<f>i,j,k)] + " -f A^V{(I>n+i,n+\,n+i) (2.4) 

where A is the cubic lattice spacing, and the conjugate momenta 
Lattice sites are located at i,j,k = 1 ,..., -f- 1 in a total sample volume V = = 

N^A^. Our phenomenological equations of motion for the lattice fields are then 
given by: 

^id,k = —A~^dH/d(f>i,j,k — F<^,-,j,fe + Ci,3,k (2.5) 


and we now have < >= 0 and < Ctj,it(OC/,m,n(<') >= 2T ASii6jnSkJ{t ~ t')/A^. 

If the noise terms are drawn from a Gaussian distribution, the Fokker-Planck 
equation [10], which describes the probability p of measuring a given configuration 
fc}, has the stationary solution: 


p oc e 


-F/A 


( 2 . 6 ) 


where F is the discrete version of the usual Ginzburg-Landau free energy functional 
(= f d^x[(V(^)^/2 -|- Vt’(^)])' Therefore, we henceforth identify the amplitude A with 
the temperature T. The damping parameter F does not explicitly appear in the <j> 
distribution at long times, which should not be surprising as the equilibrium dis- 
tribution should be independent of the dissipation time scale. The most probable 
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field configuration obtained from dFfd<j>ij^k = 0, is when all the lattice fields 

j = Tfy/r for T < Tc, and = 0 for T > 2^. A second-order phase transition is 
therefore expected near Tc. 

Before proceeding further, we find it very convenient to define, and use exclusively, 
the following dimensionless fields/paxameters: P = 4>/t),T = Tfm,i = tm^ A = Am, 
T = ttiTIti^, and R = dP/di. [In addition, we use units where h = c = ks = I 
throughout.] We solve the two first-order equations associated with Eqn. (2.5) via 
Eulerian integration. The difference equations, in our dimensionless units, are: 


(2-7) 

where G{J,k = ~ Pi+l,j,k ~ Ri-l,j,k ~ Px,j+l,k ~ Pi,j-l,k ~ Pi,j,k +1 ~ Pij,k-l]/ ^ 

is the time step, and ^ ^ ^ ^tj»/:(0^^»Tn,n(^ d" hc^ ^ — 2cXT^/,o^t7^jm^i:n/^ * 

We further define a dimensionless amplitude through the relation o:,-,j,jfe(f) = 

'^i,j,k{^) \J 24T’f c/ A®. To satisfy the discretized versions of Eqns. (2.2) and (2.3), we 
assign at each lattice site, every time step, a randomly chosen listing Q from a 
library of 50000 Gaussian distributed numbers with < Q >= 0, < >= 1/12. 

A similar Langevin approach can be found in Ref. [11], where the two-dimensional 
XY model, a model for both liquid- helium films and smectic- C liquid crystals, is 
studied. In particle physics, Langevin techniques have been used as an alternative to 
Monte Carlo simulations of lattice field theories [12]. In this case, the Langevin time 
has no physical meaning-it is simply an artificial extra dimension used to evolve field 
configurations to a state of thermal equilibrium. These techniques have also recently 
been used in cosmology to study the effects of fluctuations on the dynamics of the 
inflaton [13]. In these studies, as is the case here, we can view the Langevin time as 
‘real’ time. 

We begin by mapping out the phase structure of the field <f> on the lattice. Since 
the correlation length ^ oc l/^H (see, e.g. [14]), near the phase transition, we choose 

a temperature-dependent lattice spacing A = Nowhere shall we be concerned 

with variations in (f> on scales smaller than the coherence length so the cutoff used 
here is both computationally sound and physically reasonable. To examine phase 
structure, we define an order parameter tp- 


\ 


rP=\{P)\ 


(2.9) 
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which describes a completely broken state if V’ = 1> and the state of full [4> ^ —<p) 
symmetry i{ ip = 0. Although <p is indeed an order parameter, it is more convienent 
to use Ip as it provides a vacuum-independent description of the transition (e.g., 


depending on the initial conditions of 
negative). 

Thermal averages of local variables (' 
and translational invariance: 



run, < $ > could be either positive or 
5. P, R) are obtained from time averages, 

( 2 . 10 ) 

( 2 . 11 ) 


Measurements should be taken on time scales imax » f In addition, because 
information travels at a finite speed c = 1 in this model, we further require tmax » 
NA so that all regions in the sample volume will have communicated with each 
other for a sufficiently long time (this condition is not essential, however, if the initial 
conditions are judicously chosen) . We have taken e = A/6 and f = 3y/r, throughout. 
We have also chosen the case Tc = 0.1, which is equivalent to choosing A = 1.5 x 10“^. 
A 20 X 20 X 20 lattice, with periodic boundary conditions, was used to explore ip and 

< <p"^ > as a. function of temperature. To quickly establish equilibrium, we initially 
chose equal to the expected mean field value, and took Ri,j,k = 0. [The initial 
conditions are irrelevant in determining the long-time statistical properties, but other 
initial conditions were used as a check.] The results are shown in Figs. 2, 3. Typically, 
equilibrium configurations were reached after • 100 time steps, but several thousand 
time steps were taken in a given run. Uncertainties, and averages, were obtained from 
statistically independent “measurements” obtained from binning the data into bins of 
width the correlation time. From Fig. 2 we see that ip appears to be a continuous order 
parameter, and as expected indicates a second-order phase transition with Te 0.1. 
The mean field theory result P = ^/T {T < Tc) is also plotted in Fig. 2, in agreement 
with our numerical results ouside of the critical region. In Fig. 3 we have plotted 

< (P— < P >y > f < P as a function of the temperature. 

The equilibrium properties of domain walls wais also examined as a function of 
temperature. Wall segments were located by checking if <pij,k changed sign between 
neighboring lattice sites. All wall segments must either connect to form closed sur- 
faces, or open surfaces terminating on the lattice boundary [we refer to a given surface 
as a domain wall). The total domain wall area in the sample volume was calculated 
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as a function of temperature, and is shown in Fig. 4. Averages, and error baxs, of 
the total wall area were calculated from three statistically independent field config- 
urations. The domain wall area rapidly increases near the transition temperature. 
The transition from the broken phase to the symmetric < (f> >= Q phase, from a 
topological viewpoint, is then due to the rapid increase of domain wall area near the 
transition temperature-eventually they saturate space and < <^ >= 0 is achieved. 
The number of domain walls, of a given size, as a function of temperature was also 
calculated. Below T = 0.07 we found no walls in our runs. For 0.09 < T < 0.094 
we found 0-13 walls, in a given run, of the minimum allowed size (6 sides). As the 
temperature wais increased beyond this (to T = 0.0995), the very small walls became 
more numerous, but at the same time the wall area in smafi walls compared to that 
in large walls became increasingly smaller. 

We have ceased describing the wall distribution beyond T = 0.0995 as it is sensible 
to describe structure, e.g. a domain wall, only when the fluctuations in the field 
<f> are not too large. We, therefore, concentrate on examining structure when the 
fluctuations first start becoming unimportant below T<, (the Ginzburg criterion [15]) 
and where the mean field theory approach starts becoming applicable. Fluctuations 
in (f> have an insignificant effect on the thermodynamics of the phase transition if [15]: 

T » r<? (2.12) 

where Tc = mTc/t)^. We refer to the temperature when the above condition is near 
saturation as the Ginzburg temperature. This temperature, which is essentially de- 
fined when the change in free energy AF between the two vacua in a correlation 
volume is comparable to the temperature, serves as a benchmark for the temperature 
Tst at which structure becomes well-defined as the fluctuation rate between vacua 
is r fi^ct oc and below the Ginzburg temperature the rate is exponentially 

supressed. A naive comparison of the expansion rate Texp of the Universe with the 
fluctuation rate, where the exponential prefactor is constructed using dimensional 
analysis, indicates that structure at, or not too far below, Ta is essentially “frozen 
in” (there is still tension in the walls, however, which will smooth out structure be- 
low the freeze-out temperature T,t, see, e.g.. Ref. [16] ). In what follows, we shall 
assume that Tat < 7g, and only explore equilibrium properties below the Ginzburg 
temperature; we mention, however, that it might be possible for the undetermined 
dimensionless prefactor to alter our conclusion. 

To be more precise about the Ginzburg temperature, which by its very nature is 
somewhat indeterminate, we shall define Tq by a condition equivalent to AF ^ T, 
yet expressed in terms of statistical properties of <^, which are directly obtainable 
from our simulations: 
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(W-W)">/(#^1 (2.13) 

The above definition is very appropriate, as the criterion (2.12) can be derived by 
comparing the width of the probability distribution p{<t>corr), where <f>corr are values 
of <j> obtained in a correlation volume, with the average value of 4>. Since our lattice 
spacing is in fact a correlation length. Fig. 3 can be used to determine the Ginzburg 
temperature; Tq Ci 0.0994. We performed three runs on a 30 x 30 x 30 lattice at 
this temperature. [From Fig. 4 we note that at Tq ^ 0.0994 the total wall area is 
comparable to the maximum allowed value on the lattice, i.e., the width of the walls 
are comparable to their separation.] The distribution of the Pi,j,kS after a run is shown 
in Fig. 5. It can be seen that the width of the distribution is indeed comparable to 
the average value of P. The total wall area, averaged over three runs, associated with 
each wall size is shown in Fig. 6. It is clear from the figure that most of the wall 
area is associated with the smaller walls. An unusual feature of this graph are tiny 
peaks, first appearing for a wall area of 16^^, surrounded by much larger ones. This 
can be explained by the fact that configurations with surface area ^ (2 + 
where V is the volume enclosed by the domain wall, first start appearing when the 
wall area is 16^^ (the configuration corresponding to a minimization of surface area 
of 4 connected cubes). 

The first simulation of domain walls, in a cosmological context, was performed by 
Harvey et al. [17]. They randomly assigned + and - signs to the cells of a lattice, 
and examined the resulting structure. An identical non-dynamical simulation of the 
formation of domain walls was also performed by Vachaspati and Vilenkin [5]. When 
comparing results, we shall compare with Ref. [5], as they made the interpretation 
that the lattice spacing should roughly be a correlation length (based upon work done 
by Kibble [2]), leading to the presently accepted scenarios of wall and string formation. 
It is well known from percolation theory [18] that if lattice sites are “marked” by an 
X with probability px, then a cluster of X's will percolate through the lattice for 
px above a critical probability pc. For a cubic lattice pc ^ 0.31, which means that 
both + and - clusters will percolate in such a simulation, and that large-scale domain 
walls will always exist. The results from such a simulation are indisputable if there 
are practically no correlations in the sign of <j> on scales « However, it is clear from 
our previous work that in thermal equilibrium, at least for temperatures T < Tq, 
either -}- or - signs, hut not both, percolate through the lattice and an interpretation 
of the lattice spacing in Ref. [5] as a correlation length at the Ginzburg temperature 
is inappropriate. 

The lack of both -|- and — percolation at the Ginzburg temperature can be seen 
from mean-field theory. In the Gaussian approximation we have: dp/dP = exp{—[P — 
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PYI2W^)I\/^W , where W'^ =< {P — Py >, and P =< P >. Choosing the vacuum 
state with P > 0, the probability of finding a + site is then: 

p+ = [1 + erf{P/V2W)]/2 (2.14) 

where erf{x) is the error function. Applying condition (2.13), we find p+ = [1 + 
er/(l/x/5)l/2 ~ 0.84, and hence p_ ~ 0.16. This indicates that the — sites will 
not percolate through the lattice at (or below) the Ginzburg temperature. It is 
important to note that a transition from small-scale domain structure to large-scale 
domain structure occurs at a temperature Tl$ not too much larger than our definition 
of the Ginzburg temperature (2.13). We estimate where this transition occurs by 
setting p+ = 1 — Pc — 0.69 in Eqn. (2.14), finding that P/W ~ 0.5. The relationship 
between P/W and temperature, or AF/T, is easily estimated by neglecting the spatial 
gradient term in the free energy and fitting a Gaussian, centered at P, to 
where 1^ is a correlation volume. The result is that P^/W'^ = Xj]*t"^V^/ZT = 8AF/T, 
where AF is the change in free energy between the P = P and P = 0 states. [To 
consider AF between the two distinct vacua, the spatial gradient term cannot be 
ignored. Estimating the gradient term in the free energy as {</)— < 4> >yV^/2^'^, we 
find P^/W^ ~ 3AF/2T.] Since AF/T w P'^/W'^, we see that at T[,s the exponential 
supression in the fluctuation rate has really not “kicked in” yet, and our definition of 
the Ginzburg temperature is an appropriate, conservative one. [Recall that it is the 
freeze-out temperature which is really relevant, and which we expect to be below our 
definition of the Ginzburg temperature.] It is clear that Tis, expressed in terms of 
parameters in the potential, must have the same scaling as the Ginzburg temperature. 
This fact has also been recently noticed based upon other considerations [19]. 

It should not be a surprise that a simulation of the type of Vachaspati and Vilenkin 
cannot describe the equihbrium distribution of domain walls, even though the only 
length scale that arises in the simulation is In general, the correlation length is not 
a measure of a length over which correlation volumes are uncorrelated with respect to 
the sign of ^it is just a measure of the typical scale of spatial fluctuations of the field. 
In addition, if one simply assigns -f-’s and -’s at random, with equal probability, the 
+ — symmetry is manifest, indicating that such a simulation can only make sense if 

T > Tc, or if the fields are not completely thermalized (e.g., if there are horizons). The 
present situation indicates a fundamental problem with the Monte Carlo simulation- 
it is too simplistic. In equilibrium, at the Ginzburg temperature, aU the structure is 
essentially small-scale. However, we also know that large-scale structure must result 
due to causality being limited to horizon scales. If the correlation length and horizon 
scales are disparate, as was originally intended, it seems very hard to rationalize their 


scenario. 
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B. Cooled Fields 

We now attempt to unify both (equilibrium) small-scale and (non-equilibrium) large- 
scale structure by dynamically simulating the phase transition, on time scales that 
allow horizons in our box. Structure can then be examined at Tq- There are a number 
of ways one could do this, but we choose the simplest possible procedure: (1) as initial 
conditions we take Pi,j,ki Ri,j,k = 0, (2) we set T = Tq, and (3) we input the desired 
horizon size to correlation length ratio r, from which the number of time steps of 
the simulation are calculated. So, the initial conditions represent a state of complete 
<t> —4> symmetry (^ = 0). By setting T = To, and running the simulation on time 

scales trun smaller than the box size L, local regions of broken syimnetry will develop 
on length scales ~ trun? hut the symmetry will be globally restored as the vacuum 
state in each local region is randomly determined by thermal fluctuations. 

We first provide an estimate of r at the Ginzburg temperature. The critical 
temperature of the phase transition is estimated to be Tc 2r] [9]. The correlation 
length is ~ rn ~^ / ~ ~ where tq is obtained from Eqn. 

(2.12). Assuming the transition occured during the radiation dominated era, the 
horizon size at Tq is do = 0.6g*^^^nipilTQ, where nipi = 1.22 x 10^® GeV (the Planck 
mass), and g, is the effective particle degrees of freedom {g, = 106.75 in the standard 
model). Assuming perturbative self-interactions (A < 0(1)) we have To — and 
hence r = da/io — O.lXgZ^^^nipi/Tc. For a GUT scale transition Tc ~ 10*® GeV and 
we see that r < 0(100) for A < 0(1), i.e., the horizon is probably not significantly 
larger than the correlation length for theories of greatest interest. However, we note 
for theories with transition temperatures Tc « TquTi ^ can be a very large number. 
Explorations of such theories can be found in Ref. [20], where strings form at the 
electroweak scale {Tc ~ lO^GeV) and r ~ 10*®A. 

If r is close to unity, the Monte- Carlo type simulation of Vachaspati and Vilenkin 
should accurately describe the domain wall (and string) distribution at To- [Note 
that freeze-out is determined here by a comparison of the fluctuation rate with the 
damping rate, which has been taken as a constant F = throughout. Again the 

Ginzburg temperature is relevant for providing a benchmark temperature at which 
freeze-out occurs.] Therefore, we consider the ca^e r » 1. Unfortunately, this 
requires that we use a lattice with three disparate length scales da, and L which 
significantly reduces the range of parameter space that can feasibly be studied. We 
therefore only consider the case r = 10. As domain walls are cosmologically disastrous 
over a wide range of parameters [1], we only qualitatively describe the results of the 
simulation and reserve a more detailed analysis for global strings, in the next section. 
We performed a number of simulations on both 40^ and 50^ lattices, and found in 
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each case that both + and — clusters percolated through the entire lattice resulting in 
one ‘gigantic’ domain wall which represented about 75 % of the total wall energy. In 
simulations with + and — signs randomly distributed on the lattice sites, about 90 % 
of the total wall area is associated with one large-scale domain wall. That we obtain 
a smaller percentage of the total area in the infinite wall (by ‘infinite’, we mean a 
wall whose spatial extent is comparable to the box size) is not too surprising as our 
equilibrium results indicated a strong preference towards smaller walls. 


III. Formation of Global Strings 

A. Equilibrated Fields 

We now consider the formation of global f/(l) strings (for a good review of cosmic 
strings, see [ 16 ]). These strings are described by a complex scalar field a and a zero- 
temperature potential with a degenerate set of vacuum states connected by global 
phase transformations a — + e'‘'<r, which we take to be V(cr) = A(|cr|^ — r/^/ 2 )^/ 3 !. It is 
known that such a field theory admits cylindrically symmetric string solutions [21,22] 
where the phase \ of Higgs field a (= /e‘^/v^) winds N times in the space of 
Higgs vacua as the fields are examined in configuration space through a rotation of 
27 T in the angular coordinate 6 . The string solutions with = — 1 , 1 are topologically 
stable, as the = 0 state cannot be reached by a continuous transformation of the 
Higgs fields; global strings with |A^| > 2 are unstable to decay into [A^l vortices of 
unit vorticity [ 23 ]. The classical field theory solution for the radial profile /(r) of a 
string of unit vorticity cannot be found analytically, but has the properties /(O) = 0 , 
/(oo) = T) and a simple variational calculation [ 23 ] yields a string width w « 1 . 3 m“^, 
where m is the Higgs mass. Again we shall consider a temperature-dependent effective 
potential of the form: V($) = — + Al$|'^/ 3 ! -t- constant [r = 1 — (T/Tc)^], 
where $ is the classical paxt of <r. The probability of a given configuration {$} is 
given by p oc where F is the usual Ginzburg-Landau free energy functional. 

The most probable field configuration is then given by = r/^r/2 for T < Tc, and 
|$|2 = 0 for T > T,. 

The formulcie of the last section are easily generalized if we use the representation 
$ = ($1 -1- i$2)/\/2, where and $2 are real fields. The zero-temperature Higgs 
mass has the same form as before, i.e. m = \/Xt]/\/ 2 , and we can define the same 
dimensionless parameters as in Sec. I with the following modifications: P, = $4/77, 
Ri = dPsjdt where s = 1 , 2 . The difference equations ( 2 . 7 ), ( 2 . 8 ) are then generalized 
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to: 




^hi,k{i)s + ^Ri,j,k{i)s ( 3 - 1 ) 

+PijAiUPLk(th + PLkith - r]) + (3.2) 


where Gij,k{i)a and have the same form as in Sec. I, but with an index s, 

where s = 1,2. 

As before, we start by mapping out phase structure. Here we define the order 
parameter tp as: 


= I ($/(i)/\/2)) I = 


(3.3) 


and again jp = 1 characterizes a state of maximum disorder, and ip = 0 describes a 
state with the global $ — ^ symmetry intact. We use the same lattice spacing as 

before, i.e. A = 1/ >/r , and again use a 20 x 20 x 20 lattice to explore phase structure. 
The results are shown in Figs. 7 and 8. 

The number of strings of a given length was also calculated at various tempera- 
tures. We defined a string through a given lattice face if the phase of the Higgs field 
changed by ±27 t in a closed path around the face. The total string length on the 20^ 
lattice is shown as a function of temperature in Fig. 9. Averages, and uncertainties, 
were calculated from three statistically independent field configurations. The total 
string length rapidly increases near the critical temperature, and the transition to the 
ordered t/> = 0 state can be viewed as proceeding through the increased production 
of strings. Inside each string the symmetry is restored, and when the separation of 
strings becomes indistingushable from their widths the symmetry is restored globally. 

The Ginzburg criterion for global strings is the same form as in Eq. (2.12), but 
again we shall remove the indefiniteness 2 issociated with the >> sign, and define Tq 
by the condition: 


< ($ - ($))•($ - ($)) > / ($•} ($) = [< > + < p2 > ^ 1 (3 4) 


which is equivalent to setting the change in free energy of a uniform vacuum state 
= j?^t/ 2 to a state with a $ = 0 fiuctuation of size ~ comparable to the 
temperature. From Fig. 8 we see that To 0.0987. From Fig. 9 it can be seen 
that the probability of a string passing through a randomly chosen lattice face at 
this temperature is cs: 0.05. This probability is about a factor of 6 smaller than that 
obtained from randomly assigning Higgs phases to the lattice sites. We performed 
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three runs on a 30^ lattice at this temperature. The total loop length, averaged over 
three runs, eissociated with each loop size is shown in Fig. 10. The total string length 
at Tq is in very small loops, and no infinite strings were found. Because of the low 
loop density (relative to the Monte Carlo simulations) and hmited lattice size, it was 
not considered worthwhile to fit data and find other statistical quantities as in Ref. 
[5]. For example, it would not be useful to try obtaining a fractal dimension of the 
strings when the largest loop in the simulation has only 28 segments. 

To conclude our discussion of the equilibrium properties of global strings, we com- 
pare our results with related works. Numerical studies of the equilibrium properties 
of strings [24,25] have been performed by solving the Nambu equations of motion 
for a ‘gas of strings’ in a box. Here the total energy is a fixed parameter of the 
simulation. At low energy densities (small compared to the only scale for the energy 
density where n is the energy per length of the string) it was found that an initial 
configuration of strings will chop itself up into a large number of small loops with a 
typical size of order the lattice cutoff. Infinite strings were not found below a critical 
density pc ~ Above pc, long strings appeared and as p is increased » pc most 
of the energy density is found in long strings. These results have also been obtained 
analytically [26] by quantizing the string and counting states in the microcanonical 
ensemble. Similar results were found here, quantified in terms of the temperature. 
A transition from small-scale strings to large-scale strings occurred not far above the 
Ginzburg temperature. In analogy with the domain wall case, where we analytically 
argued that a phase transition from small-scale to large-scale structure takes place at 
a temperature Tls oc Tq {Tq < Tls < Tc), we expect the same behaviour. 

Copeland et. al. [19] have also analytically studied the behaviour of loops and 
long strings, in the Abelian-Higgs model, through the phase transition. They found, 
in thermal equilibrium, that most of the energy was in long strings at a temperature 
T„ which they defined to coincide with structure formation. They also showed that T, 
was of order the Ginzburg temperature. However, T, was defined as the temperature 
above which the loop and infinite-string partition function diverged. At temperatures 
just below T., loops quickly (exponentially) dominated the percentage of the total 
energy. We disagree with their interpretation of T, as the temperature at which 
strings formed, but their results relate the same physics as that found here. 

Finally, we should remark that we have considered quite different degrees of free- 
dom than those studied in the aforementioned works. Also, we have examined global 
strings, which have long-range interactions, whereas local, Nielsen- Olesen [21] strings 
were studied in Ref. [19], and Nambu strings were studied in Refs. [24,25,26]. We 
note, however, that our results for global strings might be more general than they 
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appear. At the Ginzburg temperature the width of global strings are comparable to 
their seperation, which effectively limits the range of interaction at formation to a 
correlation length, i.e. the global strings behave like local strings. 

B. Cooled Fields 

Now, we follow the same procedure as in the previous section to examine both small- 
scale and large-scale structure. Here, = 0)i = k{i = 0)2 = 0 were chosen 

as initial conditions so that the phase of the Higgs field was initially undefined, ev- 
erywhere. Again, we only consider the case r = 10. We performed a number of runs 
on 40^ and 50^ lattices. The results of two such runs are shown in Figs. 11 and 
12. Most of the string length (> 50%) is found to occur in the smaller loops, rather 
than in large, infinite strings. This is in contreist to the results from the Monte-Carlo 
simulations, where it was found that 80% of the total string length was in infinite 
strings. Again, we can view our results in terms of the equilibrium configuration of 
strings shifting the distribution of energy to the smaller scales. 

A few words about the meaning of these results are in order. We have not exam- 
ined the dependence of the results upon our definition of the freeze-out temperature, 
which should be determined numerically. Also, a more appropriate choice of initial 
conditions would be a thermal distribution, which could be obtained using the proce- 
dure of the last subsection. Further improvements are listed in the conclusion. Such 
improvements are not warranted within the framework presented here, as a proper 
analysis of the full-blown problem should be done in an expanding spacetime. The 
important point is that we have determined that there is no current, realistic, model 
describing the formation of walls or strings when the thermal correlation length at 
freeze-out is significantly smaller than the horizon. 


IV. Concluding Remarks 

Using a simple formalism, we have been able to study both equilibrium and non- 
equilibrium properties of global strings and domain walls. We numerically verified 
the Kibble mechanism for the production of topological defects, and shed new light on 
the details of the formation process. We found that in equilibrium, for temperatures 
T < To 1 most of the energy was associated with defects of the smallest size, and there 
was no infinite structure. For a horizon size ten times larger than the correlation 
length, we found that at the Ginzburg temperature, in both the domain wall and 
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global string case, energy was shifted (relative to the Monte-carlo simulations) from 
the large scales to the small scales. We have presented results for only one choice of 
Tc, but other values were studied and the same qualitative behaviour was observed. 
The presently held view that 80% of the string length at formation is in long strings 
deserves the qualification that it is only true if the correlation length at the Ginzburg 
temperature is comparable to the horizon length. Otherwise, such a simulation would 
be highly questionable due to the assumption of random phases on correlation length 
scales. The degree of randomness is crucial in determining statistical properties. 
The equilibrium phase transition from large-scale to small-scale structure near the 
Ginzburg temperature illustrates just how sensitive the string distribution can be to 
an assumption about the phase distribution. 

It is interesting to note that the same arguments used to describe the formation 
of strings in the ‘standard scenario’ were actually also used to show that monopole 
production could be highly supressed in the early Universe [27]. We have found 
that neither of these conflicting scenarios can be completely correct. In Ref. [27] 
they showed monopoles could be suppressed due to thermal fluctuations, without 
taking into account horizons. Their results are similar to those we obtained for walls 
and strings when we examined equilibrium properties at the Ginzburg temperature- 
we found no ‘infinite’ walls or strings. However, one cannot claim that monopoles, 
strings, or walls can be cosmologically insignificant without a proper treatment of 
horizons. 

We mention that a major question remains unanswered. What happens to the 
string distribution at formation in the limit that the horizon size to correlation length 
ratio becomes very large (r — » oo), as in low temperature phase transitions {Tc « 
rup/)? This limit is certainly not equivalent to the case where just the equilibium 
configuration applies (no infinite strings); rather, it is the limit where there are finite- 
size horizons with (/d —* 0. If we assumed a power law in the distance scales for the 
ratio of the energy in small loops to that in long strings, an extrapolation of our results 
would indicate that small loops would quickly dominate the energy distribution of 
strings in this limit. A naive way to consider this limit might be to perform a non- 
dynamical Monte Carlo type simulation involving two scales, rather than one. For 
example, in the domain wall case, one might randomly assign -f- and — vacua to cubic 
cells representing horizon volumes. Then a simulation on a lattice with spacing equal 
to the correlation length, within each horizon cell, could be performed by laying down 
+ and — signs with probabilities and p_ depending upon the vacuum state of a 
given horizon and the desired temperature. One might be led to believe that the 
total area in the infinite walls at Ta is then a surface phenomenon while the total 
area associated with the finite walls is a volume effect. 
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Finally, we call attention to a number of improvements that could be made on 
the work presented here: (1) the expansion of the Universe can easily be taken 
into account and, in this regard, the potential can also be varied according to the 
temperature- time relationship of the Universe; (2) the details of the noise term can 
be improved, and physically motivated; (3) it would be useful to numerically find out 
when freeze-out occurs, rather than using the Ginzburg temperature as a benchmark; 
(4) a systematic study of statistical properties, such as the fractal dimension of in- 
finite strings at the Ginzburg temperature, would be of interest; (5) a careful study 
of the string distribution with different ratios of horizon size to correlation length at 
formation may lead to an acceptable extrapolation of the formation details for the 
case << 1; (6) lastly, it would be interesting to perform a similar analysis for 
local strings. A number of these items shall be incorporated into a future work [28]. 

Note Added. The two-scale Monte Carlo model mentioned above was recently 
examined for arbitrary freeze-out temperatures, and several values of r [29]. It is 
found that the contribution of large-scale walls to the total wall density ~ near 
the Ginzburg temperature, in qualitative agreement with the results found here. This 
seems to be a geometrical effect, in which case one might argue that large-scale 
strings contribute a fraction ~ of the total string density near the Ginzburg 
temperature. Such geometrical arguments might explain why large-scale structure 
decreased relatively less for walls than strings in simulations with the the same value 
of r. 
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Figure Captions 

1. The total string length in loops of a given size, as a function of loop size, given 
by the non-dynamical Monte Carlo model. This simulation was done on a 35^ 
cubic lattice with periodic boundary conditions. There are four ‘long’ strings 
in this particular simulation, which represent ^ 80% of the total string length. 

2. The order parameter ^ as a function of the dimensionless temperature T. Above 
Tc 0.1 the <f> *-* —<i> symmetry is restored. The solid line is the mean-field 
result P = \/t {T > Tc). 

3. The fluctuations < (P— < P >)^ > / < P as a function of the dimensionless 
temperature T. 

4. The total number of domain wall segments on a 20^ lattice as a function of 
temperature T. 

5. The distribution of fields P, j,jt in thermal equilibrium after a run on a 30^ lattice 
at the Ginzburg temperature. 

6. The total wall area in walls of a given size, at the Ginzburg temperature, as a 
function of wall size. 

7. The order parameter ^ as a function of the dimensionless temperature T. Above 
Tc 0.1 the global $ — > e’*'$ symmetry is restored. The solid line is the mean- 
field result P = y/r (T > Tc). 

8. The fluctuations < ($— < $ >)*($— < $ »>/ <$*><$> as a function 
of the dimensionless temperature T. 

9. The total string length on the 20^ lattice as a function of the dimensionless 
temperature T. 

10. The total string length in loops of a given length, at the Ginzburg temperature, 
as a function of string length. 

11. The total string length in loops of a given length on a 40^ lattice, at the Ginzburg 
temperature with horizons at ten times the lattice spacing, as a function of string 
length. 

12. The total string length in loops of a given length on a 50^ lattice, at the Ginzburg 
temperature with horizons at ten times the lattice spacing, as a function of string 
length. 
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